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ABSTRACT 


From  a  tracking  station  of  known  geocentric  location,  the  measurements 
of  the  azimuth  and  elevation  of  a  satellite  at  three  instants  of  time  tj,  t^.  t^ 
provide  tracking  data  sufficient  for  a  preliminary  determination  of  the  orbit 
of  the  satellite  about  the  earth.  From  such  data,  the  components  (in  a  geo- 
centered  inertial  coordinate  system)  of  the  position  vectors  f  ^  r^.  r^  (corres¬ 
ponding  to  tj,  t^.  t^)  of  the  satellite  may  be  calculated  and  in  terms  of  these 
components,  all  of  the  orbital  elements  may  be  expressed.  In  view  of  the 
given  positional  data  of  the  tracking  station  and  the  directional  data  of  the 
satellite's  position,  only  the  slant  ranges  p^,  p^,  p^  are  needed  to  complete 
the  determination  of  the  vectors  r^,  r^,  r^.  To  solve  for  the  slant  ranges, 
an  iterative  procedure  is  introduced  whereby  successive  systems  of  linear 
equations  are  solved  whose  solutions  rapidly  converge  to  pj,  p^,  p^.  Sub¬ 
sequently,  all  of  the  orbital  elements  are  easily  computed.  Fundamental 
use  is  made  of  the  formulas  of  Gibbs  (see  Reference)  which  provide  approxi¬ 
mate  values  of  the  ratios  ^23/^13'  ®12^®13’  which  denotes  the  area 
of  the  triangle  bounded  by  the  vectors  r.,  f.,  f.  -  r.  .  The  method  has  been 
tested  for  accuracy  by  applying  it  to  compute  the  orbital  parameters  of 
many  parameter  families  of  hypothetical  orbits.  On  comparing  the  com¬ 
puted  values  of  the  orbital  elements  with  the  corresponding  hypothetical 
values,  the  differences  were  found  to  be  exceedingly  small.  The  average 
machine  time  for  computing  all  elements  of  an  orbit  was  approximately 
0.  05  minute. 
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I.  DYNAMICAL  BASIS  OF  THE  METHOD 

A  new  method  is  introduced  in  the  present  paper  for  the  determi¬ 
nation  of  satellite  qrbits  (or  missile  trajectories)  about  a  spherical 
earth  by  means  of  angular  measurements  only. 

Let  s,  p  denote  the  vectors  from  a  tracking  station  at  time  t  to 
the  geocenter  and  the  orbiting  satellite,  respectively.  The  satellite 
position  vector  r  relative  to  a  geocentered  inertial  frame  is  therefore 
defined  by  the  vector  equation 


r  =  p  -  s  (1) 

Let  p.  p  denote  the  slant  range  of  the  vehicle  from  the  tracking 
station  and  the  unit  vector  in  the  direction  ofp,  respectively,  so  that 
Equation  (1)  may  be  written  in  the  form 

f  =  0  p  -  s  (2) 


Let  vector  and  scalar  functions  of  t..  i  =  1,  2.  3  be  denoted  simply 

'  -  A  -  A 

by  the  use  of  the  subscript  i.  Examples  are  s  (t.)  =  s.  ,  P  (t^)  =  . 

Exceptions  to  the  use  of  this  notation  will  be  the  time  intervals  t^  ~  tp 
tj  -  which  will  be  denoted  by  respectively. 


Expanding  by  means  of  a  Taylor's  series,  we  have  in  terms  of 


r^  and  its  derivatives 


rj  =  J  (t,  +  A,)  =  f,  t  ■  +  '24*  +  •• 


(3) 


On  substituting  the  dynamical  relation 


(4) 
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into  Equation  (3),  the  equation  becomes 

/  .a  A  2  V 


.  /  ^  r'l  /i/ 

= r  777] ' ^  •  • 


Vector  multiplying  the  members  of  (5)  by  r^  yields  the  equation 


J,  X  =  (f,  X  r,)  +  a,  X  r,) 


+  •  •  • 


Similarly,  ^ 

•  • •  2 

r^  X  r  j  =  (f  2  X  ”2  )  ^2  ^  *2  ^  ^ —  +  .  .  . 


Because  of  constancy  of  angular  momentum  (r  ^  x  r^  =  ^2  ^  ^2^' 
following  equation  results  from  a  combination  of  Equations  (6)  and  (7), 


In  view  of  Equation  (5),  the  magnitudes  r^  and  r2  (of  r  ^  and  r2) 
differ  by  terms  of  order  at  least  It  follows  that  the  right  member 

of  (8)  is  of  order  ~  ^  ^2’  order  if 

Consequently,  in  either  case,  and  are  both  small, 

the  approximation  results 
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(9) 


Let  A^j  denote  the  area  bounded  by  two  position  vectors  r^,  r^  and  that 
portion  of  the  orbit  between  r^,  Pj  ,  i,  j  =  1,  2,  3,  i  <  j,  and  let 
denote  the  area  of  the  triangle  bounded  by  the  vectors  f r^,  r^  -  r.. 
The  constancy  of  the  rate  of  description  of  area  assures  that  ~ 

■^23^^2’  ®12^^1  ^  ®23^^2*  accuracy  of  (9),  therefore 

depends  upon  how  closely  the  ratio  ^|2^^23  approximated  by  the 
ratio  ®i2^®23*  form  expressions  for  the  ratios 

^12^^23  derived  in  terms  of  orbit  eccentricity  and  the  eccentric 

anomaly  angles  associated  with  the  positions  of  the  vectors  r^,  r^,  Ty 
They  are 


IL 

23 


E2  -  Ej  -  e(sin  E2  -  sin  Ej) 

fij  *  ^2  “  ®  ^3  *  sin  1112) 


(10) 


sin  (E2  -  Ej)  -  e  (sin  E2  -  sin  Ej  ) 

Tirr(ir^''^Trpr“^'"FTsTn~Er^’'^'TTnT!^ 


(11) 


where 

2,  21/2  t  1 

A. .  =  ■  J -  IE.  -  E.  -  e(sin  E.  -  sin  E.)  I  , 

ij  2  I  J  »  J  ^  J 

2/1  2.1/2  r  1 

B.  .  =  ^  -I -  I  sin  (E.  -  E.)  -  e  (sin  E.  -  sin  E.)  I  , 

ij  2  I  J  1  J  1  J 

i.  j  =  1.  2,  3  ;  i  i  j. 
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Moreover 


2  ,,  2,1/2 

A..  -  B..  =  - 

ij  ij  2 


jE.  -  E.  -  sin(E.  -  E.)j  (12) 


a**  (1  -  e**) 
- 


2,1/2 


(E.  -  E.) 
J  1 


Suppose  that,  from  a  tracking  station,  angular  measurements  are  made 
of  the  line  of  sight  to  an  orbiting  vehicle  at  times  tj,  t^.  t^.  From  these 
measurements,  the  inertial  components  of  the  corresponding  unit  vectors 
Pr  Pz’  P3  determined.  The  corresponding  system  of  equations 


*‘i  ^  ^i  Pi  ’  ®i 

may  be  written,  in  which  the  six  vectors  p. ,  s^  (i  =  1,  2,  3)  are  known. 
To  solve  for  the  unknowns  Py  p^,  p^,  use  is  made  of  Equation  (9) 
together  with  the  equation  obtained  from  (9)  by  making  the  substitution 


/  *■!  ^2 
\  **3 


^  ^2  \ 

^3-  -^2/ 


where  ^^  =  ^^+^2*  ®*g-* 


(14) 


X  r. 


z 


3 


(15) 
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Since  the  vectors  r^,  r^t  r^  are  coplanar,  constants  Cj,  c^i  exist, 
such  that 


c,  ?,  t  +  Cj  =  0  (16) 

Vector  multiplications  of  Equation  (16),  first  by  r^  and  then  by  f^, 
yield  in  virtue  of  (9)  and  (15)  (for  the  equality  sign)  c  ^  A  j  =  c ^  • 

c.  =  -  c-  so  that  (16)  assumes  the  form 


^2  ^1  ■  ^^1  ^  ^2^  **2  ^  ^1  **3  "  ® 

Substituting  from  Equation  (13)  into  Equation  (l7)  yields  the  vector 
equation 

^2^1  ^1  '  ^^1  ^^^2  ^2  ^  ^1^3P3  ■  ^  ®1  ■  ^  ^^®2  ^®3 


The  system  of  equations  obtained  from  the  components  of  vector  Equation 
(18)  possesses  a  unique  solution  p^,  p^*  P3  only  if  the  determinant 

(p  j,  p^,  Pj)  ^  0»  and  not  all  of  the  components  of  the  vector  Sj  - 
(^1  +  A  2)  ^  vanish.  The  solutions  assume  the  forms 


(O,  P2.  P3) 

D 


^2  ^Pr  P2  ’ 


(19) 


in  which  A^d  =  ^2  ®  1  ”  ^3  ®2^^1®3’  P2’  P3^'  condition 

of  the  vanishing  of  a  for  non- vanishing  values  of  Aj,  ^2  implies  that 
the  three  positions  of  the  tracking  station  corresponding  to  the  times 
tf,  t2,  t^  are  collinear.  This  can  happen  on  the  surface  of  the  earth 
only  if  the  tracking  station  is  fixed  in  inertial  space  ^i.  e.  ,  s  (t^)  = 
s  (t2)  =  s  (tj)^*  The  only  points  where  this  could  occur  are  at  the  poles 
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Substituting  the  calculated  values  for  Pj,  P3  (from  (19))  into 

Equations  (13)  and  performing  the  indicated  arithmetic  operations 
results  in  the  numerical  values  of  the  components  of  the  vectors  r^, 
r^,  fj.  The  orbital  elements  in  turn  can  be  easily  calculated  in 
terms  of  these.  In  fact,  if  0.  denotes  the  angle  from  the  line  of 
intersection  of  the  orbit  plane  and  the  z  =  0  plane  to  the  vector  r. 
(oriented  so  that  Oj  <  02)  it  may  be  easily  shown  that 


0. 

1 


cos 


-1 


-(?.)j(?j  X  ^2)3  +  (^>2  <^1  ^  ^2^1  1 

TTT- 


— I — 2 - ^ — T 

X  r2)j  +  (rj  X  r2)2 


(i  =  1,  2,  3) 


(20) 


in  which  r.  denotes  the  magnitude  of  r.  and  a  subscript  j  (j  =  1,  2) 
following  a  vector  enclosed  in  parentheses  signifies  the  jth  component 
of  the  vector.  Moreover,  the  equation  of  the  orbit  may  be  written 
in  the  form 


—  =  A  +  B  cos  0  +  C  sin  0 
r 


(21) 


in  which  A  =  1/a  (1  -  e  );  B  =  e  cos  0  /a(l-e  );C  =  e  sin  0  /a  (1  -  e  ). 

P  p 

Substituting  into  Equation  (21)  the  values  for  0^,  O2,  given  by  (20)  and 
the  values  for  r^,  r2,  r^  ,  three  linear  equations  in  A,  B,  C  result 
which  can  be  readily  solved.  Moreover,  if  fl  denotes  the  angle  from  the 
x-axis  to  the  line  of  novies,  the  elements  a,  e,  0^,  n  are  expressible 
by  means  of  the  formulas 


^ - 2 - 2 - T~ 

pr  -  (B^  +  C^) 


A 

T 


e  = 


(B^  + 


0  =  tan 

P 


-1  C 


3  ' 


n  =  tan 


-1 

(fj  x  ^2)2 


(22) 
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II.  HIGHER  ORDER  APPROXIMATIONS 

The  accuracies  of  Equations  (9)  and  (15)  are  generally  improved  by 
reducing  the  time  intervals  between  successive  observations.  Equation 
(9)  becomes  exact  for  circular  orbits  if  ^  For  near  circular 

orbits  it  is  therefore  advantageous  to  choose  =  ^2* 

The  method  of  this  paper  can  be  modified  so  that  it  can  be  satis¬ 
factorily  applied  to  cases  in  which  the  time  intervals  A  j  ,  A^  are  much 
larger  than  formerly  admissable,  by  making  use  of  equations  derivable 
from  the  formulas  of  Gibbs  (see  Reference)  for  the  ratios  of  the 
triangular  areas.  In  the  terminology  of  this  paper  these  equations  are  to 
third  order  terms 


B 

B 


il 

13 


^  ^1  ^^2  ^  ^  U  (“3  ’  “^V^l  ^  ^ 

6r2  yrj 


(24) 


Equations  (9)  and  (15)  are  first  order  approximations  of  the  exact 
equations 


ri  X  r2 


*‘2  ^  ^3 


^1  *’3 


12 


3 


23 


13 


(25) 
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in  which  the  principal  parts  of  the  ratios 

employed.  It  follows  that  if  good  estimates  of  r^,  r^  can  be  found, 

the  substitutions  of  the  ratios  ®j2^®23'  ^23^^13  (24)) 

for  ^  in  the  right  members  of  Equations  (19)  will  lead 

to  much  more  accurate  positional  determinations.  The  following 
iterative  process  serves  to  determine  good  ''initial  estimates"  of 

3  3  4 

Tj,  r2.  Since  r^  -  r^  +  r^  it  follows  that  l/r^  -  \lr^  =  3rjAj/rj, 
Consequently,  the  third  terms  of  the  right  members  of  Equation 

]((  3{c  9tc 

(24)  are  of  third  order  in  the  A's.  Let  B  ^  2,  B^^*  ®23  defined 
such  that 


B 


23 

13 


B 


12 


^3 


'1 


1  ,  ^3  24  (^^^3) 


-IS - “2T  , 

Bi3  ^3  1 


1  + 


M 

6  r. 


^  (Aj  +  ^3  ) 


(26) 


Select  an  estimate  ^  of  the  slant  range  ^od  calculate  the 
corresponding  value  by  means  of  the  equation 


(1)  _ 


Pz  '  ^2 


(27) 


(1 )  *  /  *  ( 
Substitute  r2  =  r^  '  in  Equation  (26)  and  calculate  the  ratios  (^23'^  13' 

The  first  and  second  equations  o^  (19)  with  ^.  /  A  and 

^2/^3  replaced  by  (®j2^®23'  '  ^^23^^13'  ’  respectively,  yield 

computed  values  for  Pj  and  p^,  say  P|^^^  and  Thus 


( 


A.  A 
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(2)  _  Pz*  P3^ 


(2) 


(Pl* 

- D - 


(28) 


B 


13, 


in  which  D  =  (pj,  P2>  P3)  and 


-d)  - 

^  ■  ®  1  '  I  “* —  I 

Vzi  / 


♦  \d)  /r*  \(n 

X.  *2 
s,  +  I  — J- 


B 


23 


These  slant  range  values  are  then  used  to  calculate  the  corresponding 
values  means  of  the  equations 


(2)  .  I  .  (2)  - 


(2)_|.  (2)- 


P,  P,  -  .  >'2  =  P2  P2-  ^2 


(29) 


(2)  (2) 

The  values  r  r_'  '  are  then  substituted  in  Equation  (24)  to 

^  ^  (2)  (2) 

determine  ratios  which  are  denoted  by  (®j2^®23'  '  ^^23^^13'  ‘ 

This  iteration  process  is  continued  making  use  of  Equations  (28), 

(29),  (24)  until  the  differences  (®i2^®23^^^^  ”  ^^12^^23^^^  ^ 

(B23/ Bj  3)^"^^  -  (®23^®13^^^  ^  become  negligible.  The  values  of  the 

ratios  ®j2^®23’  ^23^^13’  determined  are  used  for  the  final 

determinations  of  the  slant  ranges  Pj,  Pj  and  position  vectors 

r^  r^.  r3  by  means  of  the  equations 


Pi  - 


{a,  p^,  Pj) 

- C - P, 


“®23  ^Pp  P3^ 

^  - 15 - 


B23(Pi.  P2.  <») 
^  ^2  - - 


in  which 
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®13  -  ^12  - 

=  =  ^  »2  ^  ^  - 


and 


r.  =  p.  p.  -  s.  , 

1  X  *^1  1 


=  1.2,3 


(31) 


III.  RESULTS  OF  DOUBLE  PRECISION  COMPUTATIONS  OF  ORBITS; 

SCOPE  AND  ACCURACIES  OF  THE  DETERMINATIONS 

The  method  of  orbit  determination  was  programmed  for  double  pre¬ 
cision  computation  on  a  digital  computer.  The  means  of  testing  the  accu¬ 
racies  of  the  method  consisted  of  computing  the  orbits  (all  of  the  orbit 
parameters)  of  various  hypothetical  orbits  by  using  directional  data, 
times  of  observation,  and  known  positions  of  tracking  stations,  and  com¬ 
paring  the  computed  results  with  the  corresponding  hypothetical  values. 

A  great  multiplicity  of  hypothetical  cases  were  considered  in  which  the 
elements  8,  uu,  a,  e,  6^,  G^,  ®  ^  were  selected  from  the  following 

systems  6  =:  5°,  15°,  25°,  35° . 90°;  uu  =  4°,  12°,  20°,  28° . 68°; 

a  =  3922.  33,  4443.  94 .  16,000  nautical  miles;  e  =  0,  0.0125,  0.0150, 

.  .  .  ,  0.  85;  (Gj,  G^,  G^)  =  (6°,  6.  53°,  7.  07°),  (18°,  18.  53°,  19.  07°),  .  .  .  , 

(102°,  102.530,  103.07°) . (0°,  2°,  4°) . (0°,  5°,  6°);  G  =  -24°, 

-12°,  0°,  12°,  24°,  .  .  .  ,  60°;  a  =  lO(j-l),  j  =  1 ,  2,  3 . 13.  The  sym¬ 

bols,  a,  e  have  their  usual  connotation,  8,  denote  latitude  angle  and  the 
angle  between  the  line  of  ascending  node  and  the  meridian  plane  of  the 
tracking  station  at  the  time  of  the  initial  observation;  G^  is  the  true 
anomaly  angle  of  perigee,  G^,  G^,  G^,  are  the  true  anomaly  angles  of  the 
satellite  positions  which  correspond  to  times  tj  ,  t^.  t^  and  a  is  the  in¬ 
clination  angle  of  the  orbit  plane.  In  addition  to  the  multiplicity  of  cases 
enumerated,  similar  systems  of  cases  were  considered  which  are  near 
to  a  critical  case.  The  critical  case,  in  which  the  method  is  not  applic¬ 
able,  occurs  when  the  vectors  from  the  tracking  station  to  the  three 
positions  observed  at  tj,  t^,  t^  are  coplanar.  In  each  of  the  cases  which 
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will  be  termed  a  near-critical  case,  the  tracking  station  lies  in  the  orbit 
plane  at  the  time  of  one  of  the  three  observations.  The  accuracies  of  the 
results  pertaining  to  the  near-critical  cases  are  significant  because  many 
of  these  cases  were  found  to  require  double  precision  computations.  The 
accuracies  of  the  cases  which  are  not  near-critical  are  generally  better 
than  those  of  the  near  critical  cases.  But  even  for  the  near-critical  cases, 
the  diffe rence  s  between  the  hypothetical  values  of  the  elements  and  the 
values  computed  (from  the  observations)  are  very  small.  It  will  suffice, 
therefore,  for  the  purpo^'es  of  this  report  to  describe  the  findings  which 
pertain  to  the  near-critical  cases.  Multiple  parameter  families  of  such 
cases  were  considered  in  which  the  ranges  for  the  elements  were  as 

follows:  8  =  10°.  20° . 90°;  a  =  10°,  20° . 120°;  uu  =  sin"^  (tan  0 

cot  a);  =  sin  *  (sin  S  esc  a);  0^  =  -24°,  -12°,  0°,  ....  60°;  ®2  " 

0.  53°;  0^  =  0^  -»  0.  54°;  a  =  3922.  33,  .  .  .  ,  16,  000  nautical  miles;  e  r  0  , 

0.  0125,  0.  0150 .  20. 

The  accuracies  of  the  orbit  determinations  are,  in  general,  better 
than  those  of  the  following  typical  near-critical  cases; 

Actual  Orbit  Data  Computed  Orbit  Data 

Case  (1);  a  =  30°  a  =  29.999952° 

a  =  3922.  32999  n.  mi.  a  =  3922.  2483  n.  mi. 

e  =  0.03  e  =  0.029984383 

0  =  -24°  0  =  -24.01641° 

P  P 

Pj  =  403.  359695  n.  mi.  Pj  =  403.58686  n.  mi. 

=  405.  10132  n.  mi.  =  405.09120  n.  mi. 

=  409.  15499  n.  mi.  p^  =  409.  14477  n.  mi. 

0j  =6°  ®1  =  6.000095° 

0^  =  6. 53°  02  =  6.  530084^ 

0^  =  7.07°  0j  =  7.070072° 
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Case  (2): 


Case  (3); 


Case  (4): 


a 

= 

10° 

a 

= 

10. 000069° 

a 

= 

3922.  32999  n.  mi. 

a 

= 

3922.  2591  n.  mi. 

e 

r 

0.  03 

e 

= 

0. 029986499 

Pi 

= 

914.  40702  n.  mi. 

= 

914.  385  n.  mi. 

Pz 

= 

922.  4266  n.  mi. 

P2 

922.  4048  n.  mi. 

P3 

= 

931 . 58933  n.  mi. 

P3 

= 

931 . 5673  n.  mi. 

0 

— 

-12° 

e 

-12.  01017° 

qP 

=: 

18° 

p 

9, 

= 

18.  001388° 

= 

18.  53° 

92 

= 

18.  531377° 

®3 

= 

19. 07° 

93 

= 

19. 071367° 

a 

= 

20° 

a 

= 

20.  000279° 

a 

= 

3922.  32999  n.  mi. 

a 

= 

3922.  2658  n.  mi. 

e 

= 

0.  03 

e 

0.  029988690 

^1 

= 

1  57.  92217  n.  mi. 

oi 

= 

157.  91  886  n.  mi. 

•=2 

= 

158.  42123  n.  mi. 

P2 

= 

158.  4179  n.  mi. 

"3 

= 

1  58.  98764  n.  mi. 

P3 

= 

158.  9843  n.  mi. 

e 

— 

12° 

8 

— 

11. 991499° 

p 

9, 

= 

42° 

P 

9, 

= 

42.  000836° 

92 

= 

42. 53° 

92 

= 

42.  530828° 

= 

43. 07° 

93 

s: 

43.  070818° 

a 

120° 

a 

= 

120.  000165° 

a 

= 

3922.  32999  n.  mi. 

a 

= 

3922.  3517  n.  mi. 

e 

= 

0.  03 

e 

= 

0. 030003371 

Pi 

= 

2636.  2154  n.  mi. 

p| 

= 

2636.  2390  n.  mi. 

P2 

= 

2650.  1  334  n.  mi. 

02 

= 

2650.  1571  n.  mi. 

P3 

= 

2664.  5699  n.  mi. 

P3 

= 

2664.  5938  n.  mi. 

e 

— 

12° 

8 

— 

1  1. 998896° 

p 

»i 

= 

42° 

P 

9, 

= 

41. 999988° 

= 

42. 53° 

92 

= 

42. 529991° 

«3 

43. 07° 

93 

= 

43. 069993° 

Page  13 


Case  (5); 


Ca8e(6); 


Case  (  7) 


40° 

a 

39. 999986° 

3922.  32999  n.mi. 

a 

= 

3922.  4432  n.  mi. 

0.  01 

e 

= 

0.  010020000 

779.  80816  n.  mi. 

fi 

= 

779.  80865  n.  mi. 

806.  28246  n.  mi. 

P2 

806.  28279  n.  mi. 

833.  64203  n.  mi. 

= 

833.  64220  n.  mi. 

60° 

e 

60.  118205° 

100° 

p 

*1 

= 

99. 999972° 

100. 53° 

= 

100. 529970° 

101. 07000° 

»3 

101. 069968° 

o 

0 

a 

= 

29.  999986° 

1 6000  n.  mi. 

a 

= 

15999.  892  n.mi. 

0.  05 

e 

= 

0.  049994933 

1  1 853.  595  n.  mi. 

f’l 

= 

1 1  853.  585  n.  mi. 

1 1857.  068  n.  mi. 

f>z 

s: 

1 1  857.  057  n.  mi. 

1  1 860.  802  n.  mi. 

^3 

= 

1 1 860.  792  n.  mi. 

-20° 

e 

— 

-20. 00322° 

10° 

p 

9, 

~ 

9.  999997° 

10. 53° 

^2 

= 

10.  59997° 

0 

o 

*3 

— 

11.069997° 

30° 

a 

= 

30. 000009° 

1  6000  n.  mi. 

a 

= 

1 6001 . 358  n.  mi. 

0.  30 

e 

= 

0.  30005244 

81 1  3.  3207  n.  mi. 

= 

8113.  3793  n.  mi. 

81 26.  2608  n.  mi. 

= 

8126.  3195  n.  mi. 

81  39.  9188  n.  mi. 

P3 

= 

8139.  977  n.  mi. 

-20° 

e 

— 

-20. 003833° 

10° 

p 

= 

9.99999° 

o 

0 

*2 

= 

10. 529991° 

11. 07° 

»3 

= 

1 1. 069992° 
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Case  (8): 


Case  (9): 


Case  (10); 


a 

a 

e 

P, 


e 

1 

0. 

i 

6. 


a 

a 

e 

Pi 

Pe 

P3 

e 

p 

»2 

83 

a 

a 

e 

Pi 

Pz 

P3 

e 

p 

®i 

«2 

0. 


= 

0 

0 

a 

= 

30. 000001® 

s 

16000  n.  mi. 

a 

= 

I6OO2.  093  n.  mi. 

0.  45 

e 

S 

0. 45006593 

= 

5737.  7629  n.  mi. 

Pi 

= 

5737.  8192  n.  mi. 

= 

5751. 7785  n.  mi. 

Pz 

= 

5751. 8349  n.  mi. 

= 

5766.  6526  n.  mi. 

P3 

= 

5766.  7091  n.  mi. 

20® 

e 

-19.99675® 

10® 

p 

9, 

= 

9. 999990® 

= 

10. 53® 

82 

s 

10. 529991° 

= 

0 

0 

83 

= 

11. 069992® 

— 

30® 

a 

r 

30. 000021® 

= 

1  6000  n.  mi. 

a 

= 

16003.  720  n.  mi. 

= 

0,  60 

e 

= 

0. 60008761 

= 

3294.  5488  n.  mi. 

Pi 

3294.  6052  n.  mi. 

= 

3307.  1798  n.  mi. 

P2 

= 

3307.  2364  n.  mi. 

r 

3320.  7864  n.  mi. 

P3 

= 

3320.  8432  n.  mi. 

0 

0 

1 

6 

s 

-19. 99711® 

10® 

P 

0, 

r 

9. 999987® 

= 

10. 53® 

®2 

= 

10. 529990® 

= 

. 

0 

0 

«3 

= 

11.069994® 

— 

30® 

a 

= 

29. 999965® 

= 

1 6000  n.  mi. 

a 

= 

1  5998.  423  n.  mi. 

- 

0.  80 

e 

= 

0.  79998041 

= 

41 . 443169  n.  mi. 

Pi 

= 

41 . 442664  n.  mi. 

= 

44.  263139  n.  mi. 

P2 

s 

44.  262486  n.  mi. 

= 

65.  061448  n.  mi. 

P3 

= 

65.  060481  n.  mi. 

-20® 

e 

— 

-20.  00039® 

10® 

p 

e. 

— 

10.  000003® 

s 

10. 53® 

1 

83 

s 

10.  529995® 

= 

11.07® 

83 

= 

11.069987® 
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Caee  (11):  Cl 
a 
e 


Case  (12):  a 
a 
e 


30° 

- 16000  n.  mi. 

I.  50 

12143.  289  n.  mi. 
12169.  763  n.  mi. 
12197.  176  n.  mi. 
10. 0000° 

10. 5300° 

I I.  0700° 


a 

a 


30° 

- 1 6000  n.  mi. 

20. 0000 

351905.  91  n.  mi. 
353663.  25  n.  mi. 
355449.  83  n.  mi. 
10.000° 

10. 5300° 
11.0700° 


29. 999979° 
-15998.  314  n.  mi. 

I.  5000575 
12143.  362  n.  mi. 
12169.  836  n.  mi. 
12197.  248  n.  mi. 
9. 999996° 

10. 529995° 

II. 069995° 

29. 99990° 

-15997.  779  n.  mi. 
20.002653 
351905.  91  n.  mi. 
353663.  25  n.  mi. 
355449.  83  n.  mi. 
9. 999992° 

10. 529993° 
11.069992° 
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APPENDIX 


I 


Figure  1.  Geometry  Relating  Local  Direction  Angles  of  a  Sattelite 
Position  Vector  S.  P.  to  the  Direction  Cosines,  in 

1  i 

Inertial  Coordinates,  of  the  Unit  Vector  p,  =  S.PT/  IST^I. 

I  I  I  I 

The  vectors  (i  =  1,  2,  3)  employed  in  Equation  (19)  have  been 
referred  to  a  geocentered  inertial  coordinate  frame.  Direction  angles  of  the 
vectors  p^  will,  however,  be  observed  with  reference  to  the  horizontal 
plane  through  and  a  vertical  vector  (see  Figure  1).  Let 
denote  the  angle  formed  by  the  orthogonal  projection  of  the  unit  vector 
on  the  horizontal  plane  through  and  the  east  direction  through 
(if  0  <  Xj^  <  tt/2  the  bearing  of  S.  is  northeasterly).  Let  denote 
the  elevation  angle  of  the  line  of  sight  P^  above  the  horizontal  plane 
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at  S..  The  direction  angles  will  be  observed  (measured)  with 

respect  to  the  unit  vectors  4  Cj  originating  at  directed  in  the 

east,  north,  and  vertical  directions  respectively.  Let  t'.,  ^.denote 

the  colatitude  and  longitude  of  the  location  S.  with  respect  to  the 

inertial  frame.  As  defined  in  Section  I,  s.  =  S.  0.  Let  the  inertial 

11 

components  of  -s.  be  denoted  by  x^,  y.,  z..  It  follows  from  the 
geometry  of  Figure  1  that 

X.  =  R  sin  V.  cos  ^  y.  =  R  sin  V.  sin  z.  =  R  cos  U.  (34) 

1  11^1  1  11  1 

The  unit  vectors  T).  ,  can  be  expressed  in  terms  of  inertial 
coordinates  by  means  of  the  equations 


(35) 


in  which  the  dot  above  the  symbol  denotes  differentiation  with  respect 

to  time  t.  As  t  varies,  S.  moves  (due  to  the  earth's  angular  rotation 

t  . 

rate  ^  )  in  a  constant  latitude  circle  such  that  the  components  of  -  s. 
are  given  by 

k.  =  -  R  ^  sin  l/.  sin  6^,  y.  =  R  sin  cos  z.  =  0.  (36) 

The  inertial  components  of  the  vectors  Cj  can  now  be  readily 

calculated  by  means  of  Equations  (34)  through  (36).  The  results  are 
found  to  be 

=  -  sin  6.,  ^iZ  ~  ^i3  ~  ® 

71.,  =  -  cos  U.  cos  6.  ,  71. _=  -  sin6.  cos  I/.,  71. _  =  sin  V.  (37) 

'll  1  1  iZ  1  1  'i3  1 

C-i  =  sin  cos  ^  .,  C- 1  =  sin  V.  sin  6.,  C-  ^  =  cos  u. 

^il  1  1  i2  1  i  *i3  1 


4 
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From  the  definitions  of  the  direction  angles  X.,  #i.  it  follows  that  the 
vector  p.  is  given  by 


p.  =  (cos  /i.  cos  X.)  4.  -I-  (cos  sin  X)  rj.  +  (sin  /i.)  C. 


(38) 


The  inertial  components  of  the  vector  p^  are  found  by  substituting  from 
Equations  ( 37)  into  Equation  (38)  to  be 


p.,  =  -  cos  li.  cos  X.  sin  4 .  -  cos  u.  sin  X.  cos  V.  cos  6. 

*^il  1  i  1  1  1 

+  sin  sin  cos 

p.T=  cos  11.  cos  X.  cos  6.  -  cos  U.  sin  X.  sin  6.  cos  V. 

1  1  '^1  1  1  1 

+  sin  u-  sin  v.  sin  6., 


(39) 


p,,  =  cos  IX.  sin  >.  sin  u.  +  sin  n.  cos  v. 
*^13  '  1  11 


Equations  (39)  provides  the  means  of  calculating  the  input  data  p.j 
required  by  the  method.  Knowledge  of  the  tracking  station  locations 
(in  inertial  space)  at  the  instants  of  time  tj,  t2.  t^  furnishes  the 
values  of  ,  l/,  i  =  1,  2,  3.  The  measurements  of  the  observed 
local  angles  /i.,  X.  furnish  the  remaining  required  data.  The  angle 
from  the  x-axis  (of  the  inertial  coordinate  system)  to  the  line  of 
nodes  is  given  by  -  uu. 


REFERENCE 


Danby,  J.  M.  A.  ,  "Fundamentals  of  Celestial  Mechanics," 
The  Macmillan  Company,  New  York  (1962).  For  the 
formulas  of  Gibbs,  see  page  176,  Equations  (7,3.6) 
and  (7.  3.  7). 


